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Abstract 

A new form of governing equations is derived from Hamilton's principle of least 
action for a constrained Lagrangian, depending on conserved quantities and their 
derivatives with respect to the time-space. This form yields conservation laws both 
for non-dispersive case (Lagrangian depends only on conserved quantities) and dis- 
persive case (Lagrangian depends also on their derivatives). For non-dispersive case 
the set of conservation laws allows to rewrite the governing equations in the sym- 
metric form of Godunov-Friedrichs-Lax. The linear stability of equilibrium states for 
potential motions is also studied. In particular, the dispersion relation is obtained 
in terms of Hermitian matrices both for non-dispersive and dispersive case. Some 
new results are extended to the two-fluid non-dispersive case. 
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1 Introduction 



Hamilton's principle of least action is frequently used in conservative fluid 
mechanics [1-3]. Usually, a given Lagrangian A is submitted to constraints 
representing conservation in the time- space of coUinear vectors 

Divjfc = 0, fc = 0,...,m (1.1) 

where Div is the divergence operator in the time-space. Equation (1.1) means 
the conservation of mass, entropy, concentration, etc. 

Lagrangian appears as a function of and their derivatives. To calculate the 
variation of Hamilton's action we don't use Lagrange multipliers to take into 
account the constraints (1.1). We use the same method as Serrin in [2] where 
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the variation of the density p is expressed directly in terms of the virtual 
displacement of the medium. This approach yields an antisymmetric form for 
the governing equations 



where z is the time-space variable, "star" means the transposition, and 
is the variational derivative of A with respect to jjt 

k; = (1.3) 

Equations (1.1) - (1.3) admit particular class of solutions called potential flows 
jfe = Ofejo, «fe = const, k = l,...,m 

where (fo is a scalar function. We shall study the linear stability of constant 
solutions for potential flows. 

In Section 2, we present the variations of unknown quantities in terms of vir- 
tual displacements of the continuum. In Section 3, we obtain the governing 
system (1.2) using Hamilton's principle of least action. Conservation laws ad- 
mitted by the system (1.1) - (1.3) are obtained in section 4. For non-dispersive 
case, we obtain the equations (1.1) - (1.3) in the symmetric form of Godunov- 
Fricdrichs-Lax [4,5]. Section 5 is devoted to the linear stability of equilibrium 
states (constant solutions) for potential motions. We obtain dispersion rela- 
tions in terms of Hcrmitian matrices both for dispersive and non- dispersive 
flows and propose simple criteria of stability. In Section 6-7, we generalize our 
results for two-fluid mixtures in the non-dispersive case. Usually, the theory 
of mixtures considers two different cases of continuum media. In homogeneous 
mixtures such as binary gas mixtures, each component occupies the whole 
volume of a physical space. In heterogeneous mixtures such as a mixture of 
incompressible liquid containing gas bubbles, each component occupies a part 
of the volume of a physical space. We do not distinguish the two cases because 
they have the same form for the governing equations. We obtain a simple sta- 
bility criterion (criterion of hypcrbolicity) for small relative velocity of phases. 
In Appendix, we prove non-straightforward calculations. 

Recall that "star" denotes conjugate (or transpose) mapping or covectors (line 
vectors). For any vectors a, b we shall use the notation a*b for their sca/ar 
product (the line vector is multiplied by the column vector) and ab* for 
their tensor product (the column vector is multiplied by the line vector) . The 
product of a mapping A by a vector a is denoted by Aa., h*A means 
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covector c* defined by the rule c* = {A*h)*. The divergence of a hnear 
transformation ^4 is the covector Divvl such that, for any constant vector a, 



Div(^)a = Div(Aa) 

The identical transformation is denoted by /. For divergence and gradient 

d 

operators in the time-space we use respectively symbols Div and — , where 

oz 

z* = (t,x*), t is the time and x is the space. The gradient line (column) 
operator in the space is denoted by V (V*), and the divergence operator in 
the space by div. The elements of the matrix A are denoted by a*- where 
i means lines and j columns. If f{A) is a scalar function of A, matrix 
df 

B* = — — - is defined by the formula 
oA 

*v fdfV df 



\dA)^ da) 

The repeated latin indices mean summation. Index a — 1,2 refers to the 
parameters of components densities Pq,, velocities Uq,, etc. 



2 Variations of a continuum 



Let z = ^ j be Eulerian coordinates of a particle of a continuum and 

a volume of the physical space occupied by a fluid at time t. When t belongs 
to a finite interval [to, ti], -D(t) generates a four-dimensional domain ^ in 
the time-space. A particle is labelled by its position X in a reference space 
Dq. For example, if contains always the same particles Dq = D(to), 

and we can define the motion of a continuum as a diffeomorphism from D(to) 
into D{t) 

X = Xi(X) (2.1) 
We generalize (2.1) by defining the motion as the diffeomorsphism from the 
reference space into the time-space fl occupied by the medium in the 
following parametric form 



where Z 



't = ^(A,X) 

< 

.X = 0(A,X) 
belongs to a reference space Qq. The mappings 
'A = h{t,x.) 

< 

X = i/;(t,x) 



(2.2) 



(2.3) 
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are the inverse of (2.2). Definitions (2.2) imply the following expressions for 
the differentials dt and dx 



dt 
dx 



B 



f dX\ 



(2.4) 



where 



B 



I dg_ dg_^ 
dX ' 9X 

Formulae (2.2), (2.4) assume the form 



dx = 

Prom equation (2.5) we obtain 

dx 



+ ax 



udt + F dx 



(2.5) 



where velocity u and deformation gradient F are defined by 



u = 



dcj) ( dg^ 
'dX [dX. 



F = 



50 



d4> dg f dg' 
'dXdX [dX 



Let 



't = G{X,X,e) 



(2.6) 



(2.7) 



X = #(A, X, e) 
be a one-parameter family of virtual motions of the medium such that 

G{X,X,0) = g{X,X), $(A,X,0) = 0(A,X) 

where £ is a scalar defined in the vicinity of zero. We define Eulerian displace- 
ment C = associated with the virtual motion (2.7) 



dG a* 

r = ^(A,X,0), ^ = ^(A,X,0) 



(2.8) 



We note that ^ is naturally defined in Lagrangian coordinates. However, we 
shall suppose that C is represented in Eulerian coordinates by means of (2.3). 

Let us now consider any tensor quantity represented by / (t, x) in Eulerian 

o 

coordinates and / (A, X) in Lagrangian coordinates. Definitions (2.2), (2.3) 
involve 



/ (A,X) = f{g{X,X),4,{X,X) 



(2.9) 
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Conversely, 

/(i,x) = f|^h{t,^),^|^{t,4) 



(2.10) 



Let / (A, X, e) and / (t, x, e) be tensor quantities associated with the 

~ A 

virtual motions, such that / (A, X, e) = / (t, x, e) where A, X, t, x 
are connected by relations (2.7) satisfying / (A, X, 0) = / (A , X) or 

A 

equivalently / (i , x , 0) — / (i , x)). We then obtain 

/ (A,X,£) =/ (g(A,X,£),*(A,X,£),£) (2.11) 
Let us define Eulerian and Lagrangian variations of / 

A ~ 
A d f ~ (9 f 

Sf^-^{t,^,0) and 5f = -^{X,X,0) 
Differentiating relation (2.11) with respect to £ at £ = 0, we get 

Sf^Sf-^C (2.12) 



3 Governing equations 



Consider a four-dimensional vector jo satisfying conservation law 

Divjo = (3.1) 

Actually, (3.1) represents the mass conservation law, where jo = pv, p is the 
density and v* = (1, u*) is the four- dimensional velocity vector. Let be 
scalar quantities such as the specific entropy, the number of bubbles per unit 
mass, the mass concentration, etc., which are conserved along the trajectories. 
Consequently, if jfe = ajo, 

Divjfe = ^jo = 0, k ^ l,...,m (3.2) 

Hence j^, k = l,...,m form a set of solenoidal vectors coUinear to jo. 
Hamilton's principle needs the knowledge of Lagrangian of the medium. We 
take the Lagrangian in the form 



9z ' ' 



(3.3) 
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where are submitted to the constraints (3.1), (3.2) rewritten as 

Divjfe = 0, k = 0,...,m (3.4) 
Let us consider three examples, 
a - Gas dynamics[l-3] 

Lagrangian of the fluid is L — - p |up — e {p, rj) — pll (z), where e is the 

internal energy per unit volume, 77 = ps is the entropy per unit volume, 
s is the speciflc entropy and 11 is an external potential. Hence, in variables 
jo = pv, ji = psv , z, Lagrangian takes the form 

L = ^(^-rjo) -£(i*jo,rji)-rjon(z) = A(jo,ji,z) 

where T = (1,0,0,0). 

b - Thermo capillary fluids [6,7] 

Lagrangian of the fluid is L = -p |up — e {p , Vp, r] , Vry) — pn(z). Since 
Vp = V (l*jo) and V77 — V (l*ji), we obtain Lagrangian in the form (3.3) 

^ = a^'J^' 9^'"J 



c — One- velocity bubbly liquids [8-10] 



dt dz 

where p is now the average density of the bubbly liquid and A'^ is the number 
of identical bubbles per unit volume of the mixture. We deflne again jo = pv 

and ji = A^v. By usmg — = — v = — — and A^ = rji, we 

at oz oz l*jo 

obtain Lagrangian in the form 



.^A(j„.^.j.) 



The Hamilton principle reads: for each field of virtual displacements z 
ft — > C such that C o.iT'd its derivatives are zero on dfl, 



I 



Adfl = (3.5) 
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Since variation 5 is independent of domain and measure dQ, the variation 
of Hamilton action (3.5) with the zero boundary conditions for and its 
derivatives yields 



5 j KdVt 







(3.6) 



where -— - refers to the variational derivative with respect to jjt. In particular, 



if Lagrangian (3.3) is 



A jfc, 



we get 



(5A OK ^ 



dz 
f 



OA 



di ^ 
dz 



We have to emphasize that in (3.6). the variations S should take into 
account constraints (3.4). We use the same method as in [2, p. 145] for the 
variation of density. This method does not use Lagrange multipliers since 

A 

the constraints (3.4) are satisfied automatically. The calculation of 6 j^, 
k — 0, . . . , m is performed in two steps. First, in Appendix A we calcu- 
late Lagrangian variations 6 v, S p and S (expressions (A. 2), (A. 5) and 
(A. 6), respectively). Second, by using (2.12) we obtain in Appendix B Eulerian 

A 

variations 5 jk , k — 0, . . . , m (see (B.3)) 



A 



V dz 



(DivC)/ jfe 



dz 



Let us now define the four-dimensional covector 

5^k 



K* — 



(3.7) 



Taking into account conditions (3.4) and the fact that C, and its derivatives 
are zero on the boundary equations (3.6), (B.3) yield 



5 j KdQ 



dfl 



= / E f-Div(j,K*) 

m 



ik 



dK, 



r ™ (dKk _ fdKk\ 



dKk^ 



dz J 
]Cdn 



C dn 
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Hamilton's principle yields the governing equations in the form 



dz 



= 



(3.8) 



where are given by definition (3.7). The system (3.4), (3.8) represents 
m + d + 1 partial differential equations for m + d unknown functions u, p, 
Ofe, k — 1, . . . ,m, where d is the dimension of the Q-space. Since the matrix 



dz 



dz 



(3.9) 



is antisymmetric and all the vectors j^, A; = 0, . . . , m are coUinear, we obtain 
that ilRk io = 0. 



Consequently 



(3.10) 



k=0 



and the overdetermined system (3.4), (3.8) is compatible. 



In the case ak = const, k = 1, . . . , m and A = A jo, -;r- , the system 

V dzj 

(3.4), (3.8) can be rewritten in a simplified form 



dz 



dz 

Divjo = 

We call potential motion such solutions of (3.11) that 







(3.11) 



(3.12) 



where (fo is a scalar function. In this case, the governing system for the 
potential motion is in the form 



SA 






dz 


Divjo 


= 



(3.13) 



For the gas dynamics model, the system (3.12), (3.13) reads 



f dfo 
dt 

V(po = u* 



/ 1 d \ 



I I + div(pu) = 
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Eliminating the derivative we obtain the classical model for potential 

flows [2] 



dt 



( d\x* 



+ V 



dt 

u* = V(/?o 



d_ 
dp 



e{p) + n = 



(3.14) 



4 Conservation laws 



Equations (3.4), (3.8) can be rewritten in a divergence form. The demonstra- 

tion is performed for Lagrangian A which depends only on , — — and z. 

az 

The following result is proved in Appendix C. 



Tlicorcm 1.1 



Let L = A ( jfe , — — , z ) . The following vector relation is an identity 
az J 

Div ff: (k* j,/ - j,K* 



\k=0 



f)\ rn m 

+ 7^ + EK:Divj, -Ej^i?, ^ (4.1) 

fe=0 k=0 

dA 

where the matrices are given by (3.9) and = 



9^ 



5z 



In particular, we get 



Theorem 4.2 



The governing equations (3.4), (3.8) are equivalent to the system of conserva- 
tion laws 

Divg(Ka./-j.K:.^:|)-A/) + |^ = o 

Divjfe = 0, A; = 0, . . . ,m 
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In the special case of potential flows (3.13) the governing equations admit 
additional conservation laws 

Div (cK*o - (K*c)/^ = (4.2) 
where c is any constant vector. Indeed, 

Div(cK;) = c-(^)-. Div(K;c/, = |(K;c, = c.(^) 

Since for the potential flows 

dKp ^ /dKo\* 
dz \ dz J 

we obtain (4.2). In particular, if we take c = 1 for the gas dynamics equa- 
tions, we get conservation laws (3.14). Multiplying (4.1) by jo and taking into 
account the identity (3.10), we obtain the following theorem 



Theorem 4.3 



The following scalar relation is an algebraic identity 

(oiv (E(K*j,7 - j,K* + ^) - Ajjj jo 

BA 

+ n^Jo + E (KUo)Divjfe = (4.3) 

k=0 

This theorem is a general representation of the Gibbs identity expressing that 
the "energy equation" is a consequence of the conservation of "mass", "mo- 
mentum" and " entropy" . Examples of this identity for thermocapillary fluids 
and bubbly liquids were obtained previously in [6,10]. 

Identity (4.3) yields an important consequence. Let us recall the Godunov- 
Friedrichs-Lax method of symmetrisation of quasilinear conservation laws [4- 
5] (see also different applications and generalizations in [11-12]). We suppose 
that the system of conservation laws for n variables q has the form 

df 

-^ + divF, = 0, /, = /.(q), F, = F,(q), i = l,...,n (4.4) 

Let us also assume that (4.4) admits an additional "energy" conservation law 

Oe 

— + divE = 0, e = e(q), E = E(q) 
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which is obtained by multiplying each equation of (4.4) by some functions 
and then by summing over i — 1, . . . ,n 

1^ + divE = + divF.j (4.5) 

In particular, if we consider e, E and Fj as functions of fi, we obtain from 
(4.5) 

de _ , dE _ ^ dFi _ de dF, 

Let us introduce functions N and M such that 

TV = fip'-e, M = - E (4.7) 
Consequently, from equations (4.6) and (4.7) we find that 

dN ^ dm ^ ^ 

Hence, substituting (4.8) into (4.4), we get a symmetric system in the form 

d'^N ( g^M df \ _ ^ 

dp'-d'pj dt \ d'p'd'pj (9x ) 

If matrix Nij = ^ . „ . is positive definite then the symmetric system (4.9) 

op'op^ 

is t-hyperbolic symmetric in the sense of Friedrichs. Obviously, matrix ^ . 

is positive definite if and only if matrix e*-' = — r is positive definite, 

df'df^ 

since e^^Njp = Sp, where dp is the Kronecker symbol. 
Now, let us rewrite identity (4.3) in the non-dispersive case 

Div E (K*jfe I - jfeK*) - A / jo + E jo) div j, = (4.10) 

\fc=0 / / k=0 

Identity (4.10) is exactly of the same type as identity (4.5). It means that 
the system (3.4), (3.8) can be always rewritten in a Godunov-Friedrichs-Lax 
symmetric form. Actually, if we take 

m 

e = J2 Kfejfc - iikKDl - A and E = -Pj^K* 

fe=0 

where 

'0 10 0' 
P = I 1 
1 
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the conjugate variables p are 



— KUo 
P 



k = 0, ... ,171 



Consequently, the Gibbs identity (4.10) gives directly a set of conjugate vari- 
ables. Therefore, we have proved the following theorem 



Theorem 4.4 



If L — A(jjk), the system (3.4), (3.8) is always symmetrisahle. 



The property of convexity, needed for the hyperbolicity of the governing sys- 
tem, should be verified for each particular case. For example, in the gas dy- 
namics, the energy of the system is given by the formula 

e = E K*j, - E (j.K^)! - A = eM + 



Since 



fc=0 



KoJo 



fe=0 



u 



p 



9£ 

dp 



^ dr] 



the conjugate variables are (see also [4]) 

/ u \ 



de |u| 
'd~p 2~ 

de 

\ df] / 



P 



V 



T 



I 



where p is the Gibbs potential and T is the temperature. Obviously, if e(p, t]) 
is convex, the total energy e is convex with respect to pu, p and r]. 



5 Stability of equilibrium states for potential flows. 



5jfc 

We assume that L is a function of jfe and --— besides it does not depend 

az 

on z, i.e. 
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Let us give some definitions. An equilibrium state is a solution of equations 
(3.4), (3.8) such that = j^e = const. Let u he a real unit vector. Any 
vector f3 can be represented in the form 

^ ^ oju + (3^, where i/* = 

The equihbrium state }ke 7^ is linearly stable in the direction u if and 
only if all non-trivial solutions of the form e^l^ ^ of the system (3.4), (3.8) 
linearized at the equilibrium state jjte are such that uj is real for any real /3(T. 



5.1 Non-dispersive case 



We note that in a non-dispersive case the stability means hyperbolicity of 
governing system [13-14]. We omit the index "0" and rewrite system (3.13) in 
the form 

d] dz (5.1) 
Divj = 

where A = A(j). The Legendre transformation of A(j) is 

A(K) = K*j - A(j) (5.2) 

If the matrix A"(j) = ^ I f ^1 I is non-degenerate, (5.2) involves the 

J- = ^ (5.3) 



formula 



Hence, relations (5.1) - (5.3) yield 
~ dK [[dKl I 9^ 



A"(K)— = A"{K)^"{z) (5.4) 



where 



A"(K) = 



dK \\dKj ^ 
The second equation (5.1) and (5.4) involve 



d_ U dip 



tr(A"{K)if"{z)) = 

If we replace ip by i^e'*^ ^, we get immediately the dispersion relation 

/3*A"(Ke)/3 = (5.5) 
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where index "e" corresponds to the equihbrium state. The matrix A"(Ke) 
can be easily calculated in terms of Lagrangian A(j). Indeed, 

It follows that 

A"(Ke) = (A"(je))' (5.6) 
Relations (5.5) - (5.6) imply the following result. 



Theorem 5.1 {criterion of stability for potential non- dispersive motions) 

If the symmetric matrix G — A"(je) has the signature (— , +, +, +), the 
equilibrium state je is stable in any direction v belonging to the intersection 
of two cones 

Ci = u*G-^u < O} and C2 = {i^ | u*Gu < o} 

Proof . Since Ci and C2 contain the eigenvector corresponding to the neg- 
ative eigenvalue of G, Ci H C2 7^ 0. We can find orthogonal coordinates 
{t,xi,X2,X3) such that the dispersion relation (5.5) takes the form 

3 

-t'^ + = 0> with Aj > 

1=1 

In coordinates {t,xi,X2, xs) the cone C2 is defined by the inequality 

3 ™2 

-^^ + E r < 

i=l 



Let us represent /3 in the form 
/3 ^ 



( * ) 






( 


Xi 


= UJ 


ni 


+ 


X2 


n2 








V 



The dispersion relation /3*G /3 = imphes 

2 



/ 3 \^ 3 

\i=l ) i=l 



It involves 



-1 + E -^i^i + 2u; ^ (1 + Xi) myi + ^ ^iVi - E ^iVi 



i=l 



1=1 



a=l 
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Due to the following inequality, 



3 /3 \2 3 /3 

i=l \i=l / i=l \i=l 



2 n / -x \ 2 

\ Hi 



cu is real if u belongs simultaneously to Ci and C2. The theorem is proved. 

For the gas dynamics model, if the volume energy is a convex function of p 
and T) — ps, the matrix G satisfies the conditions of Theorem 5.1. 



5.2 Dispersive case 



If L = A I j, — ), the governing system is 
dz 



_ 5A dA 



dA 



Divl g( dl 
dz 



dip 
dz 



[Divj = 
The linearised system in a coordinate form is 



■'ifesj + i ksJ,p ^ks J 



mp 



and j' = 



(5.7) 



where the comma denotes the derivative with respect to z'' while A^^, F^^, 
D^f , calculated at point je, are defined by the relations 



A 



d^A 



d^A 



ks 



QjkQjs^ s dj^dj^p dj^dj^p 
The following symmetry relations result from (5.8) 

Aks = Ask, ^^ks = 



D 



mp 



d^A 



^ sk-> ^ks — ^sk 



(5. 



(5.9) 



For the solution of (5.7) in the form f — J^e'^"^", </? = i$e*^f^^, we get 
the following dispersion relation 



( G 



det 



/3\ 



-/3*Adj(G')/3 = 



(5.10) 



where Adj {G) is the adjoint matrix to G, and the elements Gks of the matrix 
G are defined by 



Gks — ^ks + ^r|^/3p + D]^f(3rnPp 



(5.11) 
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Equations (5.9) - (5.11) imply that G is hermitian matrix. In particular, if 
detG 7^ then Adj (G) = det(G) G^^ and the dispersion relation (5.10) is 
equivalent to (3*G~^f3 = 0, which is a generalization of (5.5) for the dispersive 
case. We get the following obvious result: 



Theorem 5.2 {criterion of stability for potential dispersive motions) 

^ks — ^'^^ symmetric matrix G defined by (5.8), (5.11) has the 
signature (— , +, +, +) for (3 = 0, the equilibrium state je is stable for small 
(3 in any direction v belonging to the intersection of the cones Ci, i = 1,2 
defined in theorem 5.1. 



We note that F^^ defined by (5.8) are always zero if the expansion of La- 
grangian A in Taylor series at the vicinity of equilibrium state does not contain 

linear terms with respect to 

az 

Sometimes, we are able to obtain a "global" stabihty. For example, let us 
consider a particular case of bubbly liquids, for N/ p = const, defined in 
Section 3 



L = -p|up + 



dp 
'di . 



e{p) 



where e{p) is a convex function of the density and a is a positive function of 
p. Then, 



L = 



- rj + 



- ^(i*j) 



Since the governing equations are invariant with respect to the Galilean trans- 
formation 

t' = t, x' = x + Ut, u' = u + U 

wc can always assume that Ug = (it is sufficient to consider the governing 
system in the reference frame moving with the velocity U = Ue). Hence, 
je = Pel, where 1* = (1,0,0,0). Omitting index "e", we can calculate the 
matrix G defined by (5.11) 



G = 



V 











I 



Hence, 



Adj(G) = 



/I 




0* 
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and dispersion relation (5.10) reads 



p3 



< - 



+ 



9p2 







(5.12) 



Let u be the direction of time in the time-space, then u = I and (5.12) is 
equivalent to 



UJ 



(5.13) 



— au 



The graph of the dispersion relation (5.13) for positive values of uo is presented 
in Figure 1. We have denoted by cu^ = \ —— - the eigenfrequency of bubbles 



Q2^ 

and by Cq = \\ P'q~^ the equilibrium sound sj^eec? of bubbly hquid (see [8]). 




Fig. 1. Dispersion relation for bubbly liquids 



6 Governing equations for mixtures 



Let us consider homogeneous binary mixtures. The mixture is described by 
the velocities Uq,, the average densities pa and the specific entropies Sa 
for each component [a = 1,2). We introduce two reference frames associated 
with each component in the form 



t — 9a{^a , Xq,) 
X = (paiK , Xa) 



(6.1) 
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and the inverse mappings 



' Ac = ha{t, x) 



(6.2) 



The corresponding famihes of virtual motions generated by (6.1) are defined 
by 



with 



_ X $q(Aq ; Xfj , £q) 

Gai^a , Xq , 0) = Qai^a , ^a) 



_ *a(A«, X„, 0) = 0„(A„, X„) 

We define the two Eulerian displacement Ca — i'^a , ^a) where 



dGr 



(Aa , Xq, , 0), 



(Aq , Xq , 0) 



As in Section 3 wc define tensor quantities associated with the two virtual 

A ~ 

motions and variations 5a and 5a- In the general case, we have two four- 
dimensional solenoidal vectors jo(a) = Pa Vq- corresponding to the a*'^ compo- 
nent, where Pa is the density and v* = (1, u*) is the four-dimensional veloc- 
ity vector. As in Section 3, we introduce additional physical quantities ak{a) 
associated with the four-dimensional vectors ]k{a) — ctfe(a)jo(a)) o; = 1,2, 
k{a) — 1, . . . , m{a) submitted to the constraints 



Divjfe(a) = 

For Hamilton's action in the form 



(6.3) 



a = A (jfc(a) , 



9z 



z 



the Hamilton principle reads: for each field of virtual displacements z e 
VL — > C,a such that Ca ^'^^ derivatives are zero on dfl, 



5a j AdQ 







Here 5a a is the derivative of the Hamilton action with respect to Sa, and Ca 
are the virtual displacements expressed in Eulerian coordinates z by means 
of equations (6.2). The method developed in Section 3 yields the equations of 
motion in the form 



m(a) 

k{a)=0 



\ dz \ dz 



0, K*(,) 



5A 



Si 



k{a) 



(6.4) 
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Divjfe(Q,) = 0, a — 1,2 and k{a) — 0,...,m{a) 

1 for tVifi rasp A — A ( i7,/_\ 



As in Theorem 4.1, for the case A = A ( jfc(a) , — , z j we can obtain the 



following identity 

2 m(a) 



<9jfc{a) ^ _ ^ J 



\a=l fe(a)=0 V 

Qjy 2 m(a:) 2 m{a) 

+ ^ + E E Div -EE Jfc(a) = 

a=l k{a)=0 a=l k{a)=0 



where 



1^ and = ^-(^ 



9 



dz J 

Hence, the governing equations (6.4) admit the conservation laws 



(2 m(a) / 
E E (^*k{a)iHa)I - jfc(a)Kfc(^) + 
a=l k{a)=0 \ 



^ dz J j ^ dz 

Divjfe(o!) = 0, a — 1,2 and A;(q;) = 0,...,m{a) 

We notice that for a one- velocity model, the number of scalar conservation laws 
is m + d + 1 where d is the dimension of the time-space. The number of 
unknown variables is m + d. Due to Theorem 4.3, the m + d + 1 conservation 
laws are connected by the " Gibbs identity" . In the case of mixtures, we obtain 
'^m{a) + d + 1 conservation laws for '^m{a) + 2d unknown variables. In 

a a 

the general case, the classical approach based on conservation laws, does not 
allow to obtain Rankine-Hugoniot conditions for this system. Nevertheless, in 
[15-17] we have obtained the jump conditions. The Hamilton principle provides 
these conditions without any ambiguity. 



7 Linear stability of mixtures. 



We consider the Lagrangian in the form [15-19] 
1 1 

L ^ -pi\ui\'^ + -P2\u2\'^ - W {pi , P2 , 771 , 772 , w) (7.1) 

where w = U2 — ui and rja is the entropy per unit volume of the a^^ com- 
ponent. A generalisation of (7.1) for thermocapillary fluids was also proposed 
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in [20]. If W is an isotropic function of w, Lagrangian (7.1) can be rewritten 
as follows 



A 



2.t^i Vl*j 



|jO(a)l 



J0(a) 



1* jo(a) 



W (l*io{l) , 1* jo(2) , 1* jl(l) , 1* jl(2) , 



where 



^ = 2 



J0(2) 



J0(1) 



(7.2) 



(7.3) 



l*jo(2) l*jo(l) 

Here we shall also restrict our study to the case Sa — const. To avoid double 
indices, we will denote jo(i), jo(2) by ji, j2. Hence, (7.2) and (7.3) involve 



A(ji, 



^ a=l \ ^ Ja / 



where W{pi,p2,i^) = W{pi, p2, piSi, P2S2, n). 

Omitting the tilde symbol, we consider the potential motions 

dA d(pa 
dia dz 
. Divja = 0, a = 1,2 



K* 



Let us consider the Legendre transformation of A(ji, 



(7.4) 



(7.5) 



A(Ki , K2) = ^ j;K„ - A 



a=l 



If the matrices 

_d_ n dpc 

are non-degenerate, relation (7.6) involves 



A{a/3) 



\ 



(/3a) 



i = 



Substituting (7.8) in Divja = 0, we get 



tr 



tr 



d f dA 



dz \ dKr 



(7.6) 



(7.7) 



(7.8) 



^ , ^ _d_( , 

^[^,dK,[dK^) dz 



with 



9 



dA ^ 



tr I J2 ^(7«)^7 
17=1 



= 



(7.9) 
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Hence, the equations of potential motions are 



(7.10) 



17=1 



In the equilibrium state ^(-ya) = const. Substituting cp^ by i^^e^^ ^ in 
(7.10), we obtain the dispersion relation in the symmetric form 

det /3*A(,i)/3\ _ 

''''^ l/3*A(i2)/3 /3*A(22)/3j 

which is the generalisation of the dispersion relation (5.5). However, to cal- 
culate A(^,„) in terms of A(^q) defined by (7.7), (7.9)), we have to solve the 
following system of matrix equations 



7=1 



IS. 



a/3-, Oa/3 



1, a — P 
0, ay^ P 



A simpler method is to consider the linearised system generated from (7.5) 



7=1 



dz 



= 0, Divj, = 



where A(^q;) are taken at the equilibrium state jke- Substituting j-y by J-y e*^ ^ 
and (fia by i^a we obtain the dispersion relation in the symmetric form 



det 



^A(n) A(2i) (3 0\ 

A(12) A(22) 13 

(3* 0* 

V 0* /3* 0/ 



(7.11) 



Equation (7.11) is presented in terms of the matrices A(^a) calculated directly 
from the Lagrangian. Let us consider Lagrangian (7.4). Suppose that the ve- 
locities of each component are the same at equilibrium (uie = U2e). Due to 
the invariance of the governing equations with respect to the Galilean trans- 
formation we assume, without loss of the generality, that Uie = U2e = 0. 
Suppressing the index "e" to avoid double indices, we get 



A 



(11) 



Pi 



dW 



J -11* 



Pi 



dpi 



ir 



A(2i) = A 



1 dW 



(12) 



A(22) — ( P2 — 



P1P2 djj, 



I - ir) 



dw \ i-ir 

dp ) p\ 



dpidpr 

ir 



ir 



(7.12) 
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All matrices A(„^) are diagonal. Now we are able to calculate the dispersion 
relation (7.11) which is the determinant of a square matrix of dimension 10. 
Let us denote 

dW d^W d^W d^W 

a = = Wi2 = , W22 = (7.13) 

o/jl dpi dpidp2 dpi 

Suppose that 

a > 0, wii > 0, W22 > 0, W11W22 — > (7-14) 

Taking into account (7.12) - (7.13) we obtain by straightforward calculations 
of the determinant (7.11) 

(a{pi+p2) + PiPTjPl - {a{p\wii + 2pip2«^i2+ ^2^22) + P1P2 {P2W22 + PiWiifj 

x(/52 ^ ^2 ^ ^2^ ^2 ^ (^2 ^ ^2 ^ ;32)2(^^^^^^ _ ^2^)p2^2 ^ q (^^^g) 

Let V — \, f3 — uju + A = Ta~\ - The dispersion relation (7.15) takes 

\Pa I 

the form 

(a(pi + p2)+PlP2)A^ - (a(p^Wii + 2piP2Wl2+P2^22) + PiP2(P2W22+PiWii))A^ 

+ {w22Wii-wl2)pipl = (7.16) 
The following result is proved in [21]: 



Theorem 7.1 



If potential W{pi, p2~ i-i) satisfies conditions (7.14)) o,^^ the roots A of the 
polynome (7.16) are real. 

dW 

Conditions (7.14) mean that internal energy U = W — w— — , where w = 

ow 

|w|, is convex [15-16, 21]. Hence, if the internal energy is a convex function 
and the relative velocity w is small enough, the equilibrium state is stable in 
time direction in the time-space. 

It is interesting to note that the Lagrangian of a two-fluid bubbly liquid with 
incompressible liquid phase (heterogeneous case) has the same form as (7.1). 
Indeed, in [18,19] the potential W for a bubbly liquid is proposed in the form 

W = p2£2o(p2o) - ^ m(c) |w]2 (7.17) 

where P2 = c P20, c is the volume concentration of gas phase and the index " 0" 
means the real density of gas and liquid phase. Concentration c is expressed 
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4 o 

by c — -n R N where it! is the average radius of the bubbles, N is the 

number of bubbles per unit volume of the mixture. The internal energy per 
unit mass £20 of the gas phase and the virtual mass coefficient m are known 
functions of P20 and c. By introducing the average density of the liquid phase 

Pi = Pio(l — c) with pio = const 

we can rewrite the potential (7.17) as 

w / P2P10 \ Pio fPw - Pi \ I 12 

Vpio - Pi/ 2 V pio / 

The increase of the volume fraction c producing the interaction between gas 
bubbles changes not only the coefficient m(c) but also the interfacial energy 
£i„(c) so that potential W can be generalized to the form 

W = p2£2o(p2o) - ^m(c) I w|' + Sinic) (7.18) 

Quantity ej„(c) produces an additional pressure term due to the interfacial 
effect [22]. Therefore, Lagrangian (7.1) describes not only binary molecular 
mixtures but also partial cases of suspensions. Case (7.17) is actually degen- 
erated: W11W22 — ^ 0. The introduction of the interfacial term in (7.18) 

involves the convexity condition W11W22 — wfo > provided > 0. 



A APPENDIX 



For any function / (A, X) (see definitions (2.9)-(2.10)) and its image f{t, X) 
in the {t , X)-coordinates, we obtain the following relations 

dl(dg_Y'^dl dl _ dg_d2(dg_\-' ^ dl 

dx[dX dt' dX dX dx[dX dX ^ ' ^ 



A.l Variation of the velocity 



The definition (2.6) of the velocity u yields 

d(f) (dg\ 

-dx U 



23 



In Lagrangian coordinates, perturbation of u is represented by the formula 

(9$ ( dG \'~^ 

u(A,X,£) = ^(A,X,5) (A,X,£)j 

Taking the derivative of u with respect to £ at £ = and using (A.l), we 
get 

~ _ dl_ (dg_y^ _ o ( dg_y^ _ ^ _ _ ^ ^ 

^ ~ ~dX [dXJ ~^'dx[dx) ~ di ~ ^di ~ dt ~ "^dt 

d d 

where — = — + u* v* is the material derivative. We obtain then 
dt dt ^ 

^ V = (/ - vr)^v {A.2) 
where v* = (l,u*), / is the identity tensor and 1* = (1,0,0,0). 

A.2 Variation of the deformation gradient 



In Lagrangian coordinates (A, X), the perturbation of the deformation gradi- 
ent is given by (2.6) 



F (A,X,£) - — - — — y — j 

Combination of the derivative of F with respect to £ at £ = and relation 
(A.l) gives 

^ ^ 1 o 1 2 O 

dX dX dX [dX J dX dX [dX J dX dX [dX J dX 

di Or f di dT\ ^ , , , 



A. 3 Variation of the density 

The mass conservation is represented by the formula 

pdetF = po(X) 
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The perturbation of p is in the form 

p(A,X,£) det F(A,X,£) = po(X) 

and consequently 

5 p det F + p5 (det F) = {AA) 
By using the Euler-Jacobi identity, 

5 (det F) = det Ftr{F-^ 6 F) 
relations (A. 3) and (A.4) lead to 

Sp ^ - tr [p{I - vr) (^.5) 



A.Jf. Variation of specific quantities 



Along each trajectory the specific quantities are constant. Consequently, 
Ofc = aofc(X). The perturbation of is such that (A, X, s) = aofc(X). 
Hence, 

5afc = (A6) 



B APPENDIX 



Formulae (A. 2) and (A. 5) yield directly the variation of jo 

* j. = {tr ((vr - /) I) / - (vr - /)|) j„ 

= - (M« C) /) jo (B.l) 

Let us consider the variations 5 j^, k = 1, . . . ,m with j^. = a^jo where the 
scalar fields Ofc are such that = (X). Due to (A. 6) 5 Ofc = and 
^ jfc = Ofe 5 jo- Consequently 

S h ^ - (DivC)/)jfc, A; = 0,...,m (5.2) 
Hence, (2.12) involves 

Sjfe = - (DivC)/)j. - ^C, A; = 0,...,m (B.3) 
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C APPENDIX 



Let A, Bhe two linear transformations depending on z. Let us define the lin- 

(dB \ ( dB \ 

A -Q-J such that for any constant vector field a, tr \ A -q—J ^ = 



tr A 



dBa , , 

-- — . As a consequence, we get 
oz J 

A— j 

Indeed, 

Div(ASa) = Div (A) B a + tr (^A 



(CI) 



dBa 



dz 



dB' 



(DivA)S+ tr (A — I I a 



Let us also recall the following useful formulae for any vector fields b(z), c(z) 



Div(bc*) = c*Divb + 



./ dcY 



Div(b*c/) = l-(b-c) = b'l^ + c'^ 
OZ oz Oz 

By using formulae (C.l) - (C.3), we have the following identities 

Div(K^.7) ^ K^^+j^^ 



Div(j,K*) = -K*Divjfe - 



dz 



Div U 



dz 



Div(A/) ^ -y: 

k=0 



d_/dh 



^Bi.iAl)^ + tr^Al^^^^^ 

/ \\ 

OA d fdik 



(C.2) 
(C.3) 

(C.4) 
(C.5) 



dAdh 
djk dz 



tr 



\ 



V 



d 



dik\ dz\dz 



A;=0 



dz 
dz \ dz 



dA 
'di 



By adding (C.4) - (C.7) we obtain 



Div { E (k^, / - j,K* + Al ^-^ 



, fc=0 



dA 

dz 



A A 9A 



(C.7) 



+ EK^Divj, -E Ja 



fc=0 



fc=0 



dz 



V dz 



= 



which proves Theorem 4.1. 
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